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Abstract 

We study the Abraham model for N charges interacting with the Maxwell field. On the scale 
of the charge diameter, R^, the charges are a distance e~^R^ apart and have a velocity ^/ec 
with e a small dimensionless parameter. We follow the motion of the charges over times of 
the order e~^/^i?^/c and prove that on this time scale their motion is well approximated 
by the Darwin Lagrangian. The mass is renormalized. The interaction is dominated by the 
instantaneous Coulomb forces, which are of the order e^. The magnetic fields and first order 
retardation generate the Darwin correction of the order e^. Radiation damping would be of 
the order e'^/^. 



1 Introduction 



Classical charges interact through Coulomb forces, as one learns in every course on electromag- 
netism. Presumably the best realization in nature is a strongly ionized gas, for which the Darwin 
correction to the Coulomb forces is of importance, since under standard conditions the velocities 
cannot be considered small as compared to the velocity of light, cf. 0, §65]. Thus, given charges, 
with positions Tq,, velocities Ua, charges e^, and masses rua, a = 1, . . . , N, their motion is governed 
by the Lagrangian 



1 1 1 



a = l a, (3=1 I " PI 

N 



+ J2 T~r^ \ra-rp\''^{uo,- [ro,-rp\){up- [rc^-rp])), (1.1) 

c denoting the velocity of light. The first term is the kinetic energy with a w^-correction of a 
strength m* depending on the precise model (m* = nia for a relativistic particle). The second 
term is the Coulomb potential, whereas the third term is the Darwin potential, which decays as 
the Coulomb potential and has a velocity dependent strength. 

On a more fundamental level, the forces between the charges are mediated through the elec- 
tromagnetic field. The instantaneous Coulomb-Darwin interaction is a derived concept only. To 
understand the emergence of such an interaction, in this paper we will investigate the coupled 
system, charges and Maxwell field, and we will prove that in a certain limit the motion of the 
charges is well approximated by the Lagrange equations for £d- 

Let us first describe how the charges are coupled to the Maxwell field. To avoid short-distance 
singularities, we assume that the charge is spread out over a distance i?<^, which physically is 
of order of the classical electron radius. Thus charge a has a charge distribution pa which for 
simplicity we take to be of the form 

Pa{x) = ea^{x), X e IR^, 

where the form factor Lp satisfies 

< G C~(1R3) , (^(x) = ^r{\x\) , ^{x) = for |x| > . (C) 

To distinguish the true solution from the approximation ( |1 . 1| ) , the position of a charge a in the 
coupled system is denoted by qa and its velocity hj v^, a = 1, . . . , N. The charges then generate 
the charge distribution p and the current j given by 

N N 

p{x, t) = pa{x - qa{t)) and j{x, t) = pa{x - qa(t))va(t), (1.2) 

a=l a=l 

which satisfy charge conservation by fiat. The Maxwell field, consisting of the electric field E and 
the magnetic field B, evolves according to 

d d 
c-^—B{x,t) = -VAE{x,t), c-^—E{x,t) = VAB{x,t)-c-^j{x,t) (1.3) 

with the constraints 

V ■ E{x, t) = p{x, t), V ■ B{x, t) = 0. (1.4) 
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The charges generate the electromagnetic field which in turn determines the forces on the charges 
through the Lorentz force equation 

d r 

■^{rrika^aVait)) = J d^x pa{x - qa{t)) E{x,t) +Va{t) A B{x,t) , t G H, (1.5) 

for a = 1,. . . ,N. Here is the bare mass of charge a and the relativistic factor 7^ = 
(1 — f^/c^) , which ensures \va\ < c. Note that there are no direct forces acting between the 
particles. Eqns. (|1.2| )- (|1.5| ) are known as Abraham model for charges. 
We define the energy function by 

N ^ 

n{E, B, q, v) = Y. ^hala + 2 y [^'(^) + ^'(^)]' (1-6) 



a=l 



with q = (gi, . . . , q^) and v = {vi, . . . , vn)- It then may be seen that the initial value problem 
corresponding to (p..2D-(P3D has a unique weak solution of finite energy and that H is conserved 
by this solution, compare with |^ for the case of a single particle. 

We assume that initially the particles are very far apart on the scale set by R^p. Thus we 
require, for a 7^ /?, that 

|ga(0)-g;3(0)| (1.7) 

with e > small. If particles would come together as close as R^, our equations of motion are not 
trustworthy anyhow. In addition, we require that the initial velocities be small compared to the 
speed of light, 

M0)\ ^ v^c. (1.8) 
Subject to these restrictions, in essence, the initial electromagnetic field is chosen such as to 



minimize the energy function Ti. from (|1.6|) , cf. Section |5.1| for precise statements and estimates. 
With these initial conditions, for the particles to travel a distance of order 6~^R^ it will take a 
time of order £:~^/^i?^/c, which will be the time scale of interest. Thus physically we consider slow 
particles that are far apart, and we want to follow their motion over long times. 

Next note that it takes a time of order e~^R^p/c for a signal to travel between the particles. 
This means that on the time scale of interest, retardation effects are small. If particles interact 
through Coulomb forces, as will have to be proved, the strength of the forces is of order e"^ since the 
distance is of order e~^Rip. Followed over a time span e~'^^'^R^/ c, this yields a change in velocity 
of order y/ec. On this basis we expect the orders of magnitude ( p..7|) and (|1.8| ) to remain valid over 
times of order e~^/^i?^/c. There is one subtle point here, however. The self-interaction of a charge 
with the fields renormalizes its mass. Thus in (|1 . 1|) the quantity cannot be the bare mass of 
the charge, the electromagnetic mass has to be added. 

In theoretical physics it is common practice to count the post-Coulombian corrections in orders 
of w/c relative to the motion through pure Coulomb forces. Thus the Darwin term is the first 
correction and of order (v/c)^. The next correction is of order [v/c)^ and accounts for damping 
through radiation. If we push the Taylor expansion in Section 3 one term further, one obtains 



a = 1, . . . , N. The physical solution has to be on the center manifold for ( |1.9| ). At the present 
level of precision it suffices to substitute the Hamiltonian dynamics to lowest order, which yields 
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dt \ dUa J 

Note that if the ratio Ca/mo, does not depend on a, then the radiation reaction vanishes and the 
system does not emit dipole radiation. The next order correction is {v/cY and of Lagrangian form. 
It is discussed in |^ and 

In general relativity, there is a huge effort to obtain corrections to the Newtonian orbits, which 
as a problem is similar to the one discussed here. The most famous example is the Hulse- Taylor 
binary pulsar, where two highly compact neutron stars of roughly solar mass revolve around each 
other with a period of 7.8 h |^. In this case (f/c) = 10"^. For gravitational systems there is only 
quadrupole radiation which is of order (f/c)^. To this order the theory agrees with the observed 
radio signals within 0,3%. In newly designed experiments one expects highly improved precision 
which will require corrections up to order {y/cY^. 



2 Main results 

We recall the initial conditions for the Abraham model (|1.2|) - (|1.5|) , where we set c = 1 throughout 
for simplicity. For the initial positions g° = gQ,(0) we require 

Cre-' <\ql-ql\<C2e-\ a ^ /?, (2.1) 

for some constants Ci, C2 > 0. For the initial velocities f ° = fQ,(0) we assume 

l^^°l<C3v^ (2.2) 

with C3 > 0. The initial fields are a sum over charge solitons, 

N N 

E{x,0) = E\x) = J2E,o^{x-ql) and B{x,0) = B\x) = J2 B,o^{x - q'J. (2.3) 

a=l a=l 

Here 

E^{x) = -V0^(x) + (v ■ \/(l)v{x))v and B^{x) = -v A V0„(a;) (2.4) 
and the Fourier transform of 0^ is given by 

Mk) = e^{k)/[e-{k-vf], (2.5) 

where it is understood that in we have to set e = Cq,. For this choice of data, the constraints 



( |1.4|) are satisfied for t = and therefore for all t. In case N = 1, the particle would travel freely, 
Qiit) = Q'l + ^i'^) t > 0, and the co- moving electromagnetic fields would maintain their form ( |2.3| ). 

In spirit, the bounds (|2.1|) and ( p.2| ) should propagate in time and the form (|2.3|) of the elec- 
tromagnetic fields, at least in approximation. On the other hand, for two particles with opposite 
charge one particular solution is the head on collision which violates the lower bound in ( |2.1| ). 



3 



Considerably more delicate are solutions where some particles reach infinity in finite time, [I^ . 
Thus we simply require that for given constants C^,C* > the bound 

C,e-^< sup \qa{t)-qp{t)\<C*e-\ a^p, (2.6) 
te[o,T£-3/2] 



holds, which implicitly defines the first time, T, at which ( |2.6| ) is violated. In fact ( |2.6| ) looks like 
an uncheckable assumption. But, as to be shown, the optimal T can be computed on the basis of 
the approximation dynamics generated by the Lagrangian ( |1.1[) . 

Under the assumption ( |2.6|) the velocity bound propagates through the conservation of energy. 
We define the electrostatic energy of the charge distributions as 

^stat = Ee^Q/rf=^fc|^(fc)|V2). (2.7) 

and compute the energy (|1.6| ) for the given initial data. Then 

N 

n{0) := nit = 0) = J2 ^y^alK) + £st.t + Oie) 

a=l 

with 7(t>) = (1 — v"^) We minimize the electromagnetic field energy Hf{t) = ^ J d?x [-E^(x, t) + 
5^(2;, t)] at time t for given p and j, i.e., for given positions q{t) and velocities v{t). Using ( |2.6|) it 
may be shown that 

N 

nit) > m^^-fMt)) + £,t.t + Oie). 

Since by energy conservation 7^(0) = Hit) and since the dominant contributions £^stat cancel exactly, 
we thus will continue to have the bound |fa(t)| — C^/e. (We refer to Section in Appendix A 
for the complete argument). Therefore 

sup Kit)\<C^^/^ (2.8) 

te[0,Te-3/2] 

with some constant > 0. 

As a next step we solve the inhomogeneous Maxwell equations for the fields and insert them into 
the Lorentz force equations. According to the retarded part of the fields, retarded positions ^^(s), 
s e [0,t], will show up. To control the Taylor expansion of qait) — qais) and thus of the retarded 
force, including the Darwin term, we will need bounds not only on positions and velocities, but 
also on Va and Va- Implicitly they use that the true fields remain close to the fields of the form 
evaluated at current positions and velocities. 



Lemma 2.1 Let the initial data for the Abraham model satisfy 1\ ), ( ^■^) , and (j^T^. Moreover, 
assume 

C,e-^< sup \qait) - qpit)l a^^p, (2.9) 

te[0,Te-3/2] 

for some T > 0. Then there exist constants C*, C„ > such that ^2.(\ ) and (2^) hold. In particular, 

SUPjg|-Q^2"£-3/2] I'Va 

it) \ <v < 1 for some v. In addition, we find C > and e > such that 

sup \vait)\<C6^ and sup \vait)\ < Ce'^^^ (2.10) 

tG[0,r£-3/2] tG[0,T£-3/2] 

in case that |eQ,| < e, a = 1, . . . , A^. In the estimates ( ^.dj ), ( \2. Sj ), and ( ^.IQ ), C and e do depend 
only on T and the bounds for the initial data, but not on e. 
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The proof of this lemma is rather technical and will be given in Appendix A. Using the bounds 
of Lemma we expand the Lorentz force up to an error of order e:^/^, cf. Lemma p.5| , which is 
the order of radiation damping (the Coulomb force is order and radiation damping a relative 
order e^/^ smaller). The terms up to order then can be collected in the form of the Darwin 
Lagrangian ( |1.1|) . We set 

4 2 J * 16 2 

rria = my,a + -e^^e and m„ = mb„ + — e^m^ 
6 15 

2 

with the electromagnetic mass me = | / d?k \ f{k)\ k'"^ and the Darwin Lagrangian 



1 £ ]^ iv 



,,/3=l 



for r = (ri, . . . , r^q) and -u = (mi, . . . , uat). The comparison dynamics is then 

d fdCj^\ _ dC 
dt \ dury J dr, 



It conserves the energy 



«B(r- S) ^ g + + - E_ . (2.12) 

Because of the Coulomb singularity, in general the solutions to (p.ll|) will exist only locally in time, 
the only exception being when all charges have the same sign, in which case energy conservation 
yields global existence. In the corresponding gravitational problem, for a set of positive phase 
space measure, mass can be transported to infinity in a finite time, [|10|- We do not know whether 
this can happen also for the Coulomb problem. 
We set 

g°=£-V° and = v^<, a = l,...,Ar, (2.13) 

with r° 7^ rjj for a ^ (3. Then ( p.l| ) and ( |2.2| ) are satisfied. During the initial time slip of order £~^ 
the fields build up the forces between particles and adjust to their motion. Thus during that period 
the dynamics of the particles is not well approximated by the Darwin Lagrangian and we correct the 
initial data of the comparison dynamics to the true positions and velocities only at the end of the 
initial time slip. To take into account that the comparison dynamics will have no global solutions 
in time, in general, we define r g]0, oo] to be the first time when either lim^^^- \ra{t) — Tp{t)\ = 
for some a ^ f3 or lim^^^- |rQ,(t)| = oo for some a holds for the comparison dynamics (|2.11| ). 
As our main approximation result we state 

Theorem 2.2 Let T > be fixed. Define t g]0, oo] as above and fix some 6o G]0,r[. For the 
Abraham model let the initial data be given by ( ^.ISj ) and ^■^ ). Furthermore we require \ea\ < e, 
with e = e(T, data) > from Lemma Let to = 4{R^ + C*£^^) . We adjust the initial data of the 
comparison dynamics such that qa{to) = £~^ra{£^^'^to) and Va{to) = y/£Ua{£^^'^tQ) , a = 1, . . . , N . 
Then there exists a constant C > such that for all t G [to,niin{r — 6o^T} £~^^'^] we have 

\qa{t) - £''r^{£^/H) \ < C^£, \v^{t) - y^Uaie'^H) I < C£^ a = 1, . . . , AT. (2.14) 
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Remarks (i) If we are satisfied with the precision from the pure Coulomb dynamics, then in ( p.l4| ) 
we loose one power in e. In this case, we can adjust the initial data of the comparison dynamics 
at time t = 0, and then {^J§ holds for all t e [0, min{r - 5o, T} e'^/'^]. 



(ii) In fact the initial data need not to be adjusted exactly a.t t = Iq, a. bound 
\qa{to) - e'^ra{e^^Ho)\ ~ and \va{to) - V^Ua{e^^Ho)\ ~ 
would be sufficient. 

3 Self— action and mutual interaction 

In this section we expand the Lorentz force term 

Fa{t) = I d^xpa{x-qa{t))[E{x,t) + Va{t) AB{x,t)]. (3.1) 

Since the fields {E, B) are a solution to the inhomogeneous Maxwell's equations, we may decompose 
them in the initial and the retarded fields, 



where 



E{x,t) = E^^\x,t) + E^''\x,t) and B{x,t) = B^^\x,t) + B^'\x,t), 



E^''\k,t) = cos \k\t E{k, 0)- I ^-^^^^k A B{k,0), 



M^\k,t) = cos \k\tB{k, 0) + i ^^^jf^*^ k A E{k, 0), 



\k\ 
n \k 

w 



M'Vk,t) = - f' ds cos \k\(t-s) J (k,s) + I Tci /'''*^!!! p(k,s)k, 
Jo Jo \k\ 

Jo \k\ 

cf. P, Section 4], with j{x,t) and p{x,t) from ( |1.2|) . Accordingly we can rewrite Fa{t) in as 

= J d^x paix - (x, t) + Vait) A (x, t)] 

+ J d^xp^{x-q^{t))[E^'\x,t)+v^{t) AB'^^\x,t)\ 
= Fi'\t) + Fi^\t). (3.2) 
First we consider F^\t). 
Lemma 3.1 For t G [to,Te-^/^], with to = ^{R^ + C*e-^), we have F^^\t) = 0. 

Proof: If S{t) denotes the solution group generated by the free wave equation in D^'^(IR^) © 
L^(IR^), it follows from (|2.3|) through Fourier transform that 
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where = (v9(x)f °, Vv9(x)). The analogous formula is vahd for B^^\x,t), with $^ to be 

replaced with ^^{x) = (0, A \/(p{x)). For fixed I < /3 < N and x G IR^ with |x - g{^| < t - 
assumption (C) yields [S'(t — — — fj3s)]j^(x) = for all s < by means of Kirchhoff's 

formula and Lemma [. . .j^ denoting the first component. As for t G [to, Te"^/^] and \x — qp\ > 



t — we obtain 



\x-qait)\ > |x-g°|-|g,(t)-g/3(t)|-|g^(t)-g°| >t-i?^-C*£-i-Cv^t 
> to/2-R^-C*e-^>R^ 



for e small by Lemma |2.1| , the claim follows. □ 

Turning then to Fj[\t) in (p.2| ), we write this term in Fourier transformed form and use (|1.2| ) 
to obtain 

N 

F!f\t) = e^FW(t) + Y: e„e^Fi5(t), (3.3) 



with 



Fi'-f^it) = l^ds I dk |(^(fc)|^e--^-'^"W-'^^(-)' I - cos |A:|(t - s) vp{s) + t ^'"^ k 

sin |A;|(t — s) , , , .A , 

-I — '-v^{t)^{k^vp{s))y (3.4) 

«, /? = 1, . . . , iV. The term Fj^^{t) accounts for the self-force, whereas F'^^{t) for /? 7^ a represents 
the mutual interaction force between particle a and particle (3. These both contributions are dealt 
with separately in the following two subsections. 

Before going on to this, we state an auxiliary result. 

Lemma 3.2 Let I < a, (3 < N , a ^ (3 . For t e [to,Te~^/'^] we have 
(a) - f ds J dk\0{k)\^e-'^-^'''^^'^"'^'^'^^ cos\k\{t-s)vp{s) 

dk \(p{k)\^ e~'^-^-P cos \k\T \vp - ir^k ■ vp)vp - riip] + 0{e'^''^), 

(h) I fds /rffc|y(fc)|V-^-['^-w-'^^(-)i 

Jo J \k\ 
= I dr J dk \m\"e-''<'^^ "^^k [l-^k. [rv, - \t%] - ^r^k ■ v,f] + 0(.^/^), 

(c) {-i) fds I dk |^(A;)|^e-^^-[^"(^)-^^(-)' '''' ~ v^{t) A {k A vp{s)) 
Jo J \k\ 

= i-t) dr J dk \^{k)fe-"'<-^ ^^^ji^ VaA{kA Vp) + 0{e'/^) . 

Here = Vait), etc., and ^^is = laif) — qpit). 

The proof is somewhat tedious and given in Appendix B. 
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3.1 Self— action 

For t G [to,T£-3/2] have 

■Vo,^{k^[v^-v^T])\ +0{e'"^). (3.5) 



sin 


k 


T 




k 





The rigorous proof of this relation is omitted since it is very similar to the proof of Lemma ^ 
given in Appendix B. It once more relies on the fact that we may Taylor expand 



by Lemma |2.1| , with = qa{t) etc. and t = t — s, whence 



Introducing 

^ f' ^^ M\k\^) ^-iik-v^)r^p^ J ^ I' dTcosi\k\T)e-'^''-''-^W, p e INo, 
Jo \k\ Jo 

Equ. (|3.5| ) may be rewritten as 

= hm (^-y'rfA;|^(A;)|'|t;,Jo-t;„Ji + ^(A;-t;,)t;,J2} 

+ J dk \(^{k)f |i [(1 - vl)k + (k ■ Va)Va]Io + i [{Va ■ Va)k - {k ■ Va)Va\h 

-^{k ■ - vl)k + {k ■ v^)v^]h^^ + (3.6) 

since = 0{e'^). Denote the term containing the Jp by J and the one containing the Ip by 1. 
To evaluate the limits t — oo, we can rely on the results from Section 4]. We first recall that 

j dk \(p{k)f' Jo ^ 0, j dk |<^(/c)|^ Ji —2me'Ja ^is i ^ oo, 

with 7q, = (1 — t>^) ^^"^ and rrie = ^ J dk \ip{k)fk~^. Moreover, V^Ji = —ikJ2, and therefore 

J (-2me7^)?)a + ■ V^(-2me7^) w« = -2me7^(^?)a + 7^(t;„ ■ ?)a)t;„) 

= -2m,[{l + vl)va + {va ■ Va)va) + 0{e^) (3.7) 

as t — > CX3, the latter equality according to the expansion 7^ = 1 + + 0{v^) = 1 + v"^ + 0{e^) 
and7^ = l + C(e). 

What concerns X, we know from & E that / dk \0{k)\^ klo 0, 



J dk\(p{k)f Iq ^ 2mc\va\ ""^arthlfd, J dk\(p{k)f (k ■ Va)kl2 ^ —2men( 



as t — > oo, where 



2 4 

fj,[v)z = {-^ - K'r^arth|t;|) z + (^(5i;^ - 3) + 3|w|"^ arth|w|) {v ■ z)v 

for |f I < 1 and z G IR^. Consequently, since s~^arth(s) = 1 + + + 0{s^) for s close to 
zero, it thus follows after some calculation that 

X — (iIq ■ ?)o)V^(^2me|t'a|~^arth|t'a|^ + ■ V^(^2me|t'Q|~"^arth|t'a|^ 

(2 22 2\ 14 , , 4, , , 

= ( 2 + Y^^a )"^e^^a + Y^"^e(^^a " + Cl^^ )• (3.8) 

Summarizing (|3.(j|), (^.Tp, and (^.8|), we arrive at 
Lemma 3.3 Fort E [tQ,Te~^^'^] we have 



3.2 Mutual interaction 

(r) 
1/3 



In this section we expand F^g{t) from (|3.4|) with f3 ^ a. For p G INq we have that 



A,:= j^dr j dk\m?e-''<-^^^-^^^T^ = {Atx)-' j J dxdy y^{xMy)\^^p + x - 

and 

rOO p 

Bp := dr j dk\^{k)\^e--'^-^'^^ cos{\k\T)T^ 

= {-P)W^ / dxdy(f{x)f{y)\^af3 + x -yf'"^ = {-p)Ap_i, 



as may be seen through Fourier transform. We hence obtain from Lemma 3^ that for /3 7^ a and 

tG [to,Te-3/2] 

+{va ■ v/3)Vi:Ao - v^iva " Vg)Ao + Ois'^^), (3.9) 

taking also into account that Ai = (47r)~\ thus V^Ai = 0. As a consequence of \^a/3\ = 0{e~^), 
cf. Lemma |2]l|, of assumption (C), and of Lemma |2.1| , it follows that in ( p.9| ) we have — V^Aq = 



C(£^), while all other terms are 0{e^). Since e.g. 



<Ce\ 



9 



with an obvious similar estimate for the other terms besides —V^Aq, we find from ( p.9| ) and after 
some calculation that for (5 ^ a and t G \tQ^Te~'^l'^\ 



\^al3\\ , , - / 1 \ , - 



87r|a/3| 



47r 



4vr|e 



t;/3 + 0(6 



An 



4vr|ea/3| 

2 



Finally, to deal with the lowest-order term we observe that with n = ^a/s/l^a/sl 



a/3 



4|ea^ 



dxdyip{x)ip{y)\ 



n + TTl 



— n 



(3.10) 



Defining i? = (a; — y)/|^a/3| = 0{e) for |?/| < i?^, we can expand = {n + R)/\n + i?| to 

obtain that tpi^R) = n + R — 3{n- R)n + 0{e^). As// dxdinp{x)(p{y){x — y) = 0, we hence conclude 
that the right-hand side of ( p.l0| ) is Thus we can summarize our estimates on the mutual 

interaction force as follows. 



Lemma 3.4 For (3 ^ a and t E [to, Te '^^'^] we have 



0/3 



4vr|e 



a/3 1 



8vr|ea/3| 



8vr|4a/3r 



^7r|^a/3| 



3 W - 



3(^/3 • ^g/?)' 
87r|ea/3|^ 



,a/3 



4vr|ea;3| = 



4vr|ea/3| = 



3.3 Summary of the estimates 

By (§3), and Lemma |Xg we find = Fj^-\t) for t G [to, Te-^l\ According to and 

Lemmas and ^]4| we hence have obtained the following expansion of the Lorentz force in ( |3.1|) . 
For t G [to, y^^-^^ 



Fait) 
Ga{q,v,v) 



\3 15 



16 
15 

^ 6^6/3 / ^al3 _ 1 
/3^a 



21^/3^ 



^Q/3 + 



21^/31- 



^a/3 



2|WI^ 



t'a -Vp) ■ ia^p) ^ 
?a/3 H TT 'f^/3 



le 



a/3 1 



Iea/3P 



(3.11) 



where to = 4(i?<^ + C*£~ ), ,^q,/3 = qaif) — Q/sit), Va = Va{t), and Vfs = Visit). Due to the Lorentz 
equation ^{m^alaVa) = F^it), cf. ( |1.5|) , we finally obtain the following lemma by calculating the 
right-hand side and expanding 7^. 
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Lemma 3.5 For t G [to, Te ^/^] we have 



1< a < iV, 



with Got from ^3. 1 j| j and Ma{v) the (3 x 3)-matrix Ma{v){z) = (rua + \m*^v'^)z + m* (t> ■ z)v for 



4 Proof of Theorem |2t2| 



We need to compare a solution {qa{t) , Va{t)) of (|1.2|)-(pr5D with data ( p.l3| ) to {ra{t),Ua{t)), where 
we let 

f„(t) = e-'roie^/H), u^{t) = ^u^ie^^/H), (4.1) 

and where the {ra(t),Ua(t)) are the solution to the system induced by (|2.11|) with data {r'^,u^). 
A somewhat lengthy but elementary calculation shows that {ra{t),Ua{t)) satisfy 



l<a< 



(4.2) 



cf. Lemma P]^ for the notation. Recalling that r g]0, oo] was defined to be the first time when 
either limt^r- \''^a(t) — r[^(t)\ = for some a P or limt^,— |rQ,(t)| = oo for some a holds, we 
find that ( [4.2| ) is valid for t E [0, (r — So)€~^^'^], for any 5o G]0,r[ which we consider to be fixed 
throughout. This leads to some useful estimates on the effective dynamics. 



Lemma 4.1 For suitable constants Co,C^,C > (depending on t, Sq, and the data) we have 



and 



Coe~^ < sup |f„(t) - fp{t)\ < C^e 

t&[0, (T-5o)e-3/2] 



sup |Ma(t)| < Cy/e. 

te[0, {T-5o)e-'i/'^] 



0^-1 



a 7^ /3, 



(4.3) 



(4.4) 



Proof: The bounds in (|]3|) follow from (|0| ) and the fact that |r„(t) -r^(t)| > Si and |ra(t)| < C 
on [0, r — 6q] for some 5i > 0, C > 0, by definition of r. Concerning (|4.4| ), by conservation of the 
energy Hb from ( p^.l2| ) we obtain C > 7^d(^^(0), 'u(O)) = li^j^ifit) , u{t)) > ^niaU^it) as long as the 
solution exists, in particular for t E [0,t — Sq]. □ 

To simplify the presentation, we henceforth omit the tilde and write (r, u) instead of (f, u) to 
denote the rescaled solution. Utilizing the bounds from Lemma 2J. and from ( [4.3|) , (|4.4| ), it may 
be seen after some calculation that 



N 



Got{q,v,v){t)-Go,{f,u,il){t)\ < CJ] (£=^|g^(t)-r^(t)|+55/>^(t)-M^(t)|+5|^^(t)-M^(t)|) (4.5) 

/3=1 

for 1 < a < and t G [0, Te^^/a] p ^g, (r - 6o)e-^/^] = [0, min{r - 6o, T}e~^/^]. Note that the term 
— r^l appears through comparison of ^a/j/l^a/j]^ to r^/j/lrc/jp, cf. the form of Ga in (|3.11|) . 
Next, a general (3 x 3)-matrix M{v) = a(f )id + b{v ® v) has the inverse 



M{v)-^ = a{vyHd + 



a{v)[a{v) + bv'^] 
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This remark shows \Ma{va) = 0{1) and \Ma{va) ^ — Maiua) ^| < C^\va — Ua\ for t G 
[0,mm{r - 5o,T}e-'^'\ Since |G„(g,t7,{/)| = 0{e'^) it foUows from Lemma gj, (pj), and §J) 



that 

TV 



K{t)-Uo.{t)\ < - r^(t)| +£5/2|^^(i) +0(^7/2) 
/3=1 

for 1 < a < and t G [to, min{r — Tje^^/^]. Summation over a and choosing e > sufficiently 
small this results in 

N N 

E -u„(t)| < CE (^'Ig.W -r„(t)| (4.6) 

for t G [to, min{r — 5q, T}e~^/^]. To use this basic estimate, we write da{t) = qa{t) — ra{t) as 

da{t) = da{to) + {t - to)da{to) + f {t - s)dc,{s) ds , da{t) = da{to) + f da{s)ds. 

J to J to 

We then obtain for t G [to,niin{r — 6o,T}e"^^'^] from (|46|) that 

D{t) < D{to) + {t~to)D{to) + Ce^ f\t-s)D{s)ds 

+Ce^/^ f\t - s)D{s) ds + C^, (4.7) 

Jto 

D{t) < D{to) +C6^ f D{s)ds + C6^^^ f D{s)ds + Ce\ (4.8) 

Jto Jto 



Dit) = max max Mafs)! and D(t) = max max \da(s) 

l<a<N se[to,t] l<a<N s£[to,t] 



where 

Application of Gronwall's lemma to ( |4.8| ) yields 

D{t)<C (D{to) +e^ + e^ J^' D{s) ds^ , (4.9) 
and utilizing this in (|4.7| ) implies 

D{t) < D{to) + {t- to)D{to) + C^Te + Ce-^l^{b{t^) + e^) + Ce'' f\t - s)D{s) ds. 

Jto 

Finally, (t - s) < C£-3/2 yields upon a further application of Gronwall's lemma that 

D{t) <c(D{to)+e-^/^D{to) + V^), t G [to, min{r - 5o, TK^'/']. (4.10) 



By assumption -D(to) = = -D(to). Therefore ( |4.1CI|) and ( [4.9|) imply ( p.l4| ). This completes the 
proof of Theorem |2.2| . □ 



5 Appendix A: Proof of Lemma |27l] 

This appendix concerns the proof of Lemma \2.1\ . We split the proof into three subsections. 
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5.1 Bounding the particle distances and the velocities 



We intend to use energy conservation to show (|2.8|) , and for that reason we calculate with (|27 
the field energy 



N 



N 



^T.J d'x [Elo {x - ql) + B\ (x - )] + - ^ j Sx [E,o (x - g° ) ■ (x - gj) 

a=l a,/3 = l 

,0 



+B.i{x-ql)-BAx-ql) 



According to ( p.4[ ) and Section 2] the first term equals 



a=l 



lo£ 



1 + \v). 



Denoting the term in [. . .] as %lj{r) is odd, and hence Taylor expansion implies %lj{r) 

1 + (9(r^) for r small. Therefore ( p.2|) yields 



Coul 



with 6^coui from (|2.7|). To deal with the contributions for a 7^ /3 in the second term, we obtain by 
passing to Fourier transformed form and observing ( p.2|) that e.g. 



£xE,.^{x - ql) ■ E^oix - g°) = e,e^ j d'k imfk-'e^'-^'^-'^^^ + 0{e) = 0{e), 
the latter with ( |2.1| ) and by passing to polar coordinates. Thus we have shown 

^f(O) =^Coul + 0(£). 

Next we will investigate the field energy at time t > 0. We claim that 
n^it) = lJd'x [E\x,t) + B\x,t)] >^Jd^xE\x,t) > -^{p{-,t),A-'p{-,t)] 
The easiest way to see this is to introduce potentials A and 0, 



L2(R3)- 



(5.1) 



(5.2) 



dA 

B{x,t) = V A A{x,t), E{x,t) = -V(f){x,t) ~ F{x,t), with F{x,t) = —{x,t), 

for the electromagnetic field. Then p = V ■ E = — A0 — V-F, and the estimate in ( |5.2| ) follows by 
passing to Fourier transformed form. On the other hand, substituting p from ( |1.2|) into 



l{p{.,t),A-'p{.,t))^ 



d'k\p{k,t)rk 



2l-2 



by assumption 
and thus 



we can argue exactly as before to show that the terms with a ^ (3 are 0{e) 

N 



n^t) >-J2j d'k \pMfk~' + 0{e) = fcoui + 0{6) 



(5.3) 
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for t e [0,T£-3/2], Consequently for t G [0, T£-3/2] by energy conservation, cf. ( p. .61 ), by ( |5.1j ) and 

(PD 

N N N 

E ^halK) + ^coui + 0{e) = J2 "^ba7«) + ^f(O) = J] mb„7K(t)) + ^^(t) 

> "^ba7 (^a (^) ) + ^^Coul + 0{e), 

a=l 

with -f{v) = (1 - Thus 

N N 
a=l a=l 

with some constant C depending on Ci, C3, C*, T. This estimate now allows to prove ( |2.8| ). Define 

/+ = {a G {1, . . . , iV} : ^Mt)) < 7^)} and /_ = {a G {1, . . . , iV} : ^Mt)) > ^{v'j}. 

For a G /+ we have \va{t)\ < K| < Cgx/i by dg^). Thus for e so small that C|e < 1/2, 
7«) - 7(^^a(t)) < v^KO' - (^aW)'l < Ce. Therefore by (U) 

Ce> E mb„(7(t;a(t))-7(^°))- 

Since > we deduce that 

7(^aW) <7«)+C£, aG/_, 

and according to < Csy/e it then follows that < Cy/e also for a G /-. This concludes 

the proof of 

Using (U) and 
we have 

\qa{t) - qp{t)\ < \ql - gjl + - g°| + |g^(t) - g?| < C^e-' + 2C,T^ee-"^ = C*e-\ 

with C* = C2 + 2CJT. We remark that for the estimates in this section the smallness of the Cq, 
was not needed. 

5.2 Bounding li^al^)! 
Since ^ 

with the (3 x 3)-matrices moa{va) given through moaiva){z) = mba(7a-z + lai'^a ■ z)va), z G IR^, 
we obtain from ( |1.5|) that for a = 1, . . . , iV 

Va = moaiVa)'^ J X Paix - Qa) {[E{x) - Ey^{x - Qa)] + Va A [B (x) - B^^{x - Qa)]) 

= moa{Va)~^ J d^X Pa{x)(^Zi{x + qa,t) + Va A Z2{X + qa,t)^ + Ra{t), (5.5) 
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2.8) it is finally easy to derive the upper bound in (UM), since for t G [0,Te ' ] 



where mQa{Va)~'^ z = my^a'^la^iz — {va ■ z)va), z G H^, is the matrix inverse of mQa{va)- For ( |5.5| ) 
it is important to note that adding the E^^i^x — gQ,)-term and the Va A By^{x — gQ,)-term does not 
change the integral, as may be seen through Fourier transform using ( p.4| ) and ( |2.5| ). Moreover, in 
3D we have set 



N 



13=1 



(5.6) 



and 



Z{x,t) 



Zi{x,t) 
Z2{x,t) 



E{x, t) - EjLi E^^(i^{x - qp{t)) 
Bix,t)-E^iB.^,(t)ix-qf,it)) 

Maxwell's equations and the relations (f-V)-Ei,(x) = —V AB^{x)+eip{x)v, {v-'V)Bv{x) = VAE^^x) 
e = Cq for index a, yield 

Z{t) = AZ{t) - f{t) , with A = 

and 



VA 
-VA 



(3=1 \ (Vfsit) ■ Vv)B^^{x - qp{t)) 

The Maxwell operator A generates a C^-group f/(t), t G IR, of isometries in L^(IR^)^ 
see p|, p. 435; (H2)]. Therefore we have the mild solution representation 



Z(x,t) = [f/(t)Z(-,0)](x) 



ds [U{t-s)f{-,s)]{x). 



(5.7) 
(5.8) 

L2(IR3)3; 

(5.9) 



According to ( p.3|) , Z{{)) = 0, so the first term drops out. To estimate the remaining term, we first 
state and prove some auxiliary lemmas that will be used frequently. 

Lemma 5.1 For given f = (/i,/2) with V ■ /i = and V ■ /2 = we have for W{t,s,x) = 

{Wiit, S, X), W^it, S, X)) = Pit - S)fi; S)] (X) 



Wi{t,s,x) 

W2{t,S,x) 



1 



47r(t - s)2 Jiy^^i=(^t-s) 

1 



d^y (t - s)V A /2(y, s) + My, s) + {{y - x) ■ V)/i(y, s) 



47r(t - S)2 J\y-:c\=it-s) 

Proof: See P, Lemma 8.1]. 



d'y -it- s)W A /i(y, s) + /2(2/, s) + ((y - x) ■ V)/2(y, s) 



□ 



Lemma 5.2 (a) Let ^(s) > 6e some function. Assume that for y G IR'^, s G [0,t], and some 
fiy, s) = {fi{y, s), f2{y, s)) with V ■ A = = V ■ /2 

TV ^ 

\fiiy,s)\ + \f2iy,s)\ < cas)T.r-i TT^' (5-^0) 

TV ^ 

|V/i(i/,.)| + |V/2(y,.)| < ce(5)E^ rrpj, (5-11) 

f3=il + \y-qp{s)\ 
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Then for each a = 1, . . . , N , t & [0,Te ^/^], and \x\ < 

fds[Uit-s)f{-,s)]{x + qo.it)) <C( sup as)). 
JO Vse[o,i] / 

(b) Under the hypotheses of (a), if instead of ( ]5. and i\5.1^ ) it holds for fixed 1 < a < N that 

1 



\fi{y,s)\ + \f2{y,s)\ < cas)Y: 



1 1 + b - (ll3{s)\ 



3 ' 



/3 

N 

\^fiiy,s)\ + \Vf2{y,s)\ < Cas)T., , , , „ 

l + \y-q(^{s)\ 



4' 



(5.12) 
(5.13) 



then fort G [0,Te ^/^] and |x| < we have even that 



ds [f/(t-s)/(-,s)](x + g,(t)) 



< Ci sup ^(s) )£. 

VsG[0,t] 



('cj Le^ ^(t, s) > 5e some function. Assume that for y G IR^, r E [0,t], s E [0,t], and some 
giy, 1", s) = {gi{y, r, s), 5(2(1/, r, s)) wzt/i V ■ 5-1 = = V ■ 5-2 that 



N 



\9i{y,^,s)\ + \g2iy,T,s)\ < CaT,s)J2 



1 1 + b - Qais] 



|3 ' 



\Vgiiy,T,s)\ + \\/g2iy,T,s)\ < Ce(r,s)^ 



1 1 + b - 



|4 • 



(5.14) 
(5.15) 



r/ien for each a = 1, . . . , N , t E [0,Te '^^"^j, and \x\ < -R^ 



dr / ds [U (t - s)5f(-, r, s)] (x + g„(t)) 
Jo 



<C[ sup ^(r, s; 

(T,s)eAt 



w/iere At = {(r, s) : r G [0,t],s G [0,r]}. 

/n (a)-(c), all constants C on the right-hand sides are independent of a, t, and x. 



Proof: (a) Define W as in Lemma |5.1| . We derive the estimates with Wi. Fix 1 < a < A^, 
t E [0,T£-3/2], g ^ jQ^^]^ g^^^ < According to Lemma ^ ( |CT1| ), and (§]TTp 



|iyi(t,s,x + g„(t))| <C^/^^ f:/g](t,s,a;), 

/3=1 



(t - S)2 



with 



^af]it^ Sj x) 



d'y 

\y-x-ga{t)\=(t-s) 



(t-s) 



.1 + \y - qp{s)r' I + \y - qp{s)r - 

In the sum in ( |5.16| ), with general > 2, we first consider the term lj^{t, s, x), i.e., the one with 
P = a. In this case according to ( |2.8| ), \y — q/3{s)\ > \y — x — qa{t)\ — |x| — \qa{t) — qa{s)\ > 



(5.16) 
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{t- s) - Cy/6{t - s) > {t- s)/2- for e small. Therefore \y - g^(s)| > {t - s)/4 for 

s <t — AR^. We hence obtain ioi (3 = a and s <t — AR^ 

I^-Kt.s,x)<C-^^^^^. (5.17) 
On the other hand, for s G [t — 4i?^, t] 

It2it.s,x) < Cit-sf[{t-s) + l]<Cit-sf[AR^ + l]<Cit-s)'^-^^^ 

Hence (|5.17| ) shows that the latter estimate holds for any s G [0,t]. Since 

I l + (f-s)2 i "^^L l + (t-.)3 



the term with (3 = a will satisfy the claimed estimates not only in (a), but also in (c). 

Next we turn to deriving a bound for I^p{t^ s, x) with (3 ^ a. First note that for some portion 
of the interval [0,t] the preceding argument applies again. For this, define to = 4:{R^ + C*e^^). 
Then for s < t — we find by (|2.8|) for e small that on the y-sphere 

\y-q(3{s)\ > \y - X ~ qa{t)\ - \x\ ~ \qa{t) - qi3{s)\ 

> {t-s)-R^- \q^{t) - qp{t)\ - \qp{t) - 

> {t~s)-R^- C*e-^ - Cy/e{t - s) > {t - s)/2 - R^ - C*e-^ > (t - s)/4. 
Therefore as in ( ^.17] ) for general n > 2 

/g(t,g,x)<C ^|^~'^^^„ , se[0,t-to], (5.18) 

(2) 

and it remains to estimate I^i^(t,s,x) ioi (3 ^ a and s E \t — tQ,t]. To do so, we note that an 
explicit computation shows for zi, Z2 G IR^ and 7 > 

2 1 7r7 (l + {-f+\zi-Z2\f\ 

d y — \ ?K = 1 r log 



|j/-^i|=7 {l + \y-Z2\) \Zl-Z2\ V1 + kl -^2|) 

7^7 1 /^i , 471^1-2:21 
log 1 + 



2 



2 



\zi-Z2\ \ 1 + (7 - - Zsl) 
< 2, (5-19) 

" i + (7-ki-^2|)' 



as log(l + A) < A for A > 0. Similarly, forn > 2 

J \y—zi \=-y 



-2i|=7 {l+\y-Z2\ 



2717^ 



1 dr 

-1 



1 + (\zi - Z2f + 2-f\zi - Z2\r + 72 j 
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(n+l)/2 



J-1 ^ 



dr 



Cn'i 



7 



(l + \zi - + 2^\zi - Z2\r + 
1 



(n+l)/2 



\Zl - Z2\ 



So in particular 



J\v—z^ \=-y 



2i|=7 {l+\y~Z2 



1 7 



n > 2. 



^1 - Z2 



(5.20) 



(5.21) 



Below we will also need some more refined estimates, and for this purpose we note that according 
to ( ^^201) also 



1 



'l?/-2i|=7" {I + \y - Z2 

Analogously we obtain 

1 



< C- 



1 + i\zi - Z2\ + -fY 



< c- 



I d'y 



<c 



|y-2i|=7 {l + \y-Z2\) 1 + {\zi - Z2\ - "fY 



min < 1 



Zl - Z2 



|2 • 



Zl - Z2 



(5.22) 



(5.23) 



As to bound /^^^(t, s, x) ioi (3 ^ a and s G [t— to, we then use (|5.21| ) and ( ^.19|) with zi = x+qa{t), 
Z2 = lisis), and 7 = t — s to obtain for s E [t — to, t] 



I^;^{t,s,x)<C 



+ 



X + qc^it) - g^(s)| 1 + [(t _ s) _ + g^(t) - g^(^ 



(5.24) 



Therefore by (|T8D and (IQip 



(t - s) 

< 



sup ^(s; 



ds (2), 

+ 



ds 



i!^kt,s,x) 



< C( snp m)( r\ 

\se[o,t] / \Jo 1 



ds 



ds 



+ 



+ (t - s)2 A-to |a; + g„(t) - g/3(s)| 



(5.25) 



'to 1 + [(t - s) - |x + ga(t) - g/3(s)[ 
The first of the three integrals is bounded by a constant. Concerning the second, we have 

\x + qa{t) - qp{sY > \qa{t) - q/3{tY - \x\ - \q/3{t) - qp{sY > C^e'^ - ~ Cy^{t - s) 
by ( p.6| ) and (|2.8|). In the domain of integration [t — to, t] it holds that t — s < < Ce^^, whence 
\x + qa{t)-qp{sY > C,e-^-R^-Ce-'^^ > {C,/2)e~\ s G [t-to,t], /3 ^ a, < R^, (5.26) 

for e small. Therefore the second integral can be bound by Ce Jl_^^^ ds < Ceto < C. To estimate 
the last integral =: J on the right-hand side of (|5.25|), we substitute 6 = t — s to obtain 



J 



to 



dd 



i + [e-r{eY{ 

18 



(5.27) 



with r{9) = \x + g„(t) - - 9)\. Observe that |f(^)| < - 9)\ < by dgj). Thus 

6 1-^ x(^) = d — r{6) is strictly increasing, and we can substitute 6 = 0{x) to get 



J 



x(to) dx 



1 



< C 



dx 



nl + x' 



< C. 



/x(o) l-r(e)Vl + x^ 
Summarizing these estimates we obtain the bound claimed in part (a) of the lemma, 
(b) Defining 1^^] as in (|5.16| ), we need to show 



(5.28) 



*-*o ds 



"*o ds {t — s) 



By 

In the domain of integration, {t — s) > to < Ce^'^, and hence 



Jo 1 



+ (t- sY 



< Ce. 



(5.29) 



Thus it remains to estimate the part of the integral in ( |5.28| ) for s G [t — to)^]- Firstly, by ( p.22| ). 



* ds 



-to (t - s)2 Jiy-^-g^(t)\={t-^s)'' {l + \y- g/3(s)|^) 

(j c ft 

, 1 T-^ < Ce^ \ ds = Ceho < Ce. 

-to \x + qa{t) - qp{s)\^ Jt-ta 



< c 



(5.30) 



Here we have used \x + qa{t) — ?/3(s)| > {C^/2)e ^ for e small, cf. (|5.26|) . Reference to this is 
possible, since we again have that (3 ^ a. Analogously we infer from ( ^.23| ) that 



* ds 



d'y 



(t-s) 



t-to {t - S)2 Jly-^-q^(^t)\ = {t-^s)^ (1 + |y - g/3(s)h 

(t-s) / 1 



< 



ds 

t~to \X 



+ qait) - qpis)\' VI + [{t -s)-\x + q^it) - qp{t 



< CeHn 



ds 



CeJ < Ce, 



It-to l+[{t- S) -\X + qait) - g/3(s)|]2 

with the bounded J from (|5.27|). This together with (p.30[) and ([5. 291) shows that (|5.28D is satisfied. 



(c) Due to the remarks in (a), ( |5.14| ), and ( |5.15| ) we only have to prove 



dr 



ds 



I^akt,s,x) < C, P^a, te [0,Te-^/% \x\ < R. 



(t — s 



\2 "/3 



(5.31) 



We decompose the domain of integration = {(r, s) : r G [0, t], s G [0, r]} in A^^i = fl {(r, s) : 
s G [0, t — to]} and A^ 2 = {(t, s) : r G [t — tQ,t], s E [t — to, t]}. On A^j we can utilize ( |5.18| ) to get 



drds 



I%{t,s,x) <C I dr I ds 



t,i 



(t-s) 



2 Q:/3 



1 + (t - s) = 



< C. 



(5.32) 
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Since again t-s<to< Ce^^ for (r, s) G A^^a, by (|]26D and dOTP 



< C 



- s)^ 
drds 



y-x-qo,{t)\={t-s) l + \y - qi3{s)\' 

1 1 



< Ce 



dr 



In addition, by ( |5.23|) 



At- 



< 



drds 



At,: 



1 



{t - s) \X + Qait) - qf3{s)\ 
ds rt 



t-to t 



Ce ds = CeU < C. 



(5.33) 



{t - s) J\y-x-g^(t)\={t~s) l + \y- qp{s)\' 
drds 



(t-s) 1 + [{t -s)-\x + q^it) - qp{s) 



. - - - (5-34) 
It-t, 1 + [(t _ s) - |a; + - qp{s)\Y 

since the last integral is just J from (|5.27|) and hence bounded. By ( ^.32| ), (|5.33| ), and ( |5.34|) we 
thus have proved ( |5.31| ). □ 

Lemma 5.3 Define (p^i^x) through (j)y{k) = e(p{k)/[k^ — {k-vY]. Then for x G and \v\ < v < 1, 
with V = Vx, 

|V0„(x)| + \V,VMx)\ + |V^V</>,(x)| + < C\e\{l + \x\)-\ 

I VV.(a;)| + |V,VV.(a;)| + \\/yMx)\ + |V3VV.(a;)| < C|e|(l + 

\v'u^)\ + |v.vv.(x)| + \vyM^)\ + \VlV'Mx)\ < C\e\il + \x\)-\ 
|VV.(x)| + |V,VV.(x)| + |V^VV.(x)| + |V3vV.(a;)l < C|e|(l + \x\)-\ 
Proof: Tedious calculations; see also the appendices of ^. 

Now we can estimate /q cis [[/(t - + cf. (|]9D, for t G [0,Te-^/^] and |x| < R, 

using Lemma ^T!] and Lemma |5]^(a), with / = (/i, defined by (pl8|). Since V ■ = 0, and 
V ■ -E^ = ev? is independent of v, we have V ■ /i = = V ■ /2. Concerning ( ^.10| ) and (|5.11|) , 
note \V^eJ{x)\ + \VA{x)\ < C(|V0^(x)| + |V„V0„(x)|) < C|e|(l + and |V„VE„(x)| + 

|V„VS„(x)| < C{\V^Mx)\ + |V^VV^(x)|) < C|e|(l + by Lemma p|. Thus ( CT) and 

are satisfied with ^(s) = (^maxi<^<Ar |{;/3(s)|^ (^maxi</3<Ar |e/3|). As Z{x,0) = 0, hence (|5 



□ 



in conjunction with Lemma p.2| (a) yields for a 



\Zix + qait),t)\<C{snp^me^x^\^^^^^^ t e [0,Te~'/'] 



\x\ < 



(5.35) 



We will utilize this further in (|5.5|) , and to this end we also need to bound Ra{t) from ( |5.6| ). For 
fixed j3 ^ a one calculates for the interaction terms 



d^xpaix - g«(t))V0j,^(t)(x - qfsit)) 



(— i) CaCis / d^k k 



\m\' 



Jk-[q^(t)-qa{t)] 



k^-{k-vp{t)r 



d'^xd^yipix - qa{t))f{y - qp(t))V Cvp{t){.x - y), 
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(5.36) 



with 



[{1 - v^)x^ + {x ■ vff^ 



If I < 1. 



(5.37) 



Then sup^gjQ 2-^-3/2] \^Cvp{t){x)\ < C{1 + ^ due to ( |2.8| ). By (C), in ( |5.36| ) we only need to 
integrate over {x,y) that have \x — qa(t)\ < R<fi and \y — qi3{t)\ < R^. Then by ( p.6|) , \x — y\ > 
Wait) - qf3it) \ - 2R^ > - 2R^ > {C^/2)6-^ for e small Therefore ( lOel ) shows 



By definition of By{x) and E^{x) we have 



(5.38) 



i?a(t) = moa(fa(t)) ' E ( " ^^/^W + t^/^W " ^a/?^]^/?^ + Mt) A [-^/^(t) A *a/3(t)]) (5.39) 



cf. ( p.6|) , and therefore ( ^.38| ) together with (|2.8| ) implies 

\Rait)\<C6^, t e [0,T£-3/2]. 

Hence (U), O), and (|O0| ) finally yield 



(5.40) 



K{t)\ < ^(^sup^ max^|f^(.)|)( m^ax^le^l) +6^^ 



for every a = 1,. . . ,N and t G [0,Te •^/^]. Choosing maxi<^<Ar |e^| < e with e > sufficiently 
small, we find that for a = 1 N 



sup |'ya(^)| < C^'^- 

te[0,Te-3/2] 



For later reference we also note that then according to ( ^.35|) 

\Z{x + qa{t),t)\<Ce'^, a = l,...,N, t E [0,Te-^/% \x\ < R, 

5.3 Bounding li^al^)! 
By ( p.8|) we have in particular that 

Mt)-Vfsit)\ < Cv^, t G [0,T£-3/2]. 



(5.41) 



(5.42) 



(5.43) 



In order to estimate the derivative of Equ. ( ^.5| ), first note that using the explicit form of moo(fa) ^ 
we obtain from ( ^.41|) that 



moa{Va(t)) ^ 



< C\va{t)\ < Ce\ 



Hence by (|^ and QCTI ) 



M<C{e'+\M^{t)\ + \R^{t)\), 



(5.44) 



(5.45) 



with Rn defined in (15.61), and 



M„(t)= / d=^a;p„(x) (L,(t)Zi)(x + g<,(t),t) + f«(t) A (L«(t)Z2)(x + g„(t),t) 



(5.46) 
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where La(t)(j) = {va(t) ■ V)0 + dt(p for a function = t). We first estimate Ma(t). Let 
= {Lait)Z){x,t). Since generally |[La(t)0] = Lc,(t)0+({;o-V)0 and, see (l^), ^ = AZ-f 
with / from ([5 .81), we obtain 



t^ = AJ:a + iVa-V)Z-La{t)f. 

According to ( ^.31 ) it may be shown that Tia{x, 0) = 0. We hence get 

Sa(x + g„(t), t) = ^* dr [f/(t - r) ((^,(r) • V)Z(-, r) - L,(r)/(-, r))] (a; + 

As a consequence of ^(VZ) = V{AZ — /) = ^(VZ) — V/ and 0) = 0, we obtain from the 
group property of U{-) that 

+ ga(t),t) = 



dr 



f/(t-r)((i;,(r).V)Z(-,r)) (x + g,(t)) 



— dr ds 
Jo Jo 



U{t-s){Mr)-V)f{;s))]{x + q^{t)). 



With 5f(?/, r, s) = Va{'r)-'V f{y, s) it follows from the definitions of /, E^{x), and B^{x) that V-g = 0. 
Moreover, by ( p. 411 ) and Lemma |5.3| we find that ( |5.14D and ( [5.15D are satisfied with ^{t, s) = e^. 
Therefore Lemma |5.2| (c) applies to yield for a = 1, . . . , 



\^a,iix + q^it),t)\<Ce\ te[0,Te-'/% \x\ < R^. 



(5.47) 



To estimate 



dr 



Uit-r)(L^{T)f{;r)) (x + g„(t)) 



observe that 



[L„(r)/(-,r)](x) = t;,(r)-V/(x,r) + 9i/(x,r) 



N 



J2 {(^/3 • V.)<l>^^(x - qp) + {vp ■ Vj^<l>^^(x - qp)} 
0=1 

N 



with all time arguments taken at time r, and = {E^, B^). Since V ■ -B^ = and V ■ = e^p is 
independent of we have that V ■ = = V ■ In addition, satisfies ( ^AOD and ( pi]) with 

^^(r) = ( max Ivflfr)!)! max |efl|)+e^. 
Because has an additional x-derivative, moreover ( |5.12| ) and ( |5.13| ) hold for with 



C (r) = [m.a^^\vo,{T)-v^{T)\)e 



22 



as again follows from Lemma p.3| and ( |5.41| ). Thus Lemma |5.2| (a) and (b) imply that for all 
a = l,...,N,t e [0,T£-3/2], < 



< 



dT[Uit-r)f\;r)]{x + q^it)) 



dr [[/(t-r)/^(-,r)](x + g,(t)) 



< C sup e{r)+ sup e{r)e 

Te[o,t] re[o,t] 



< c 



e^+( sup max |i;^(r)|)( m^^^ + ( sup max |t;^(r) - t;^(r)|)e^ 



Hence by (|^) and (|5^) for a = 1, . . . , A^, t G [0, Te-^l\ and |x| < K 



|S,(x + g,(t),t)| < C 



+ ( sup ^maxji)/3(r)|)( max 



According to the definition of Ma{t) in ( |5.46| ) we therefore have 



|M„(t)| 



/ ci^xpa(x)fe„,i(x + ga(t),t) +Va{t) AT,^^2{x + qa{t),t) 



< c 



e^^'^ + i sup max liiflfr)!)! max lefll 



(5.4J 



To further estimate the right-hand side of ( ^.451 ), we have to bound Ra{t), with Ra{t) from (|5.6|) . 
Calculating -RQ,(t) explicitly we obtain 

Rait) 

—moa{Vay^^moaiVa)Rait) 

^ /3=1 

+ moa(t;a)"M / d^xpa{x-qa) {{Vo, - Vp) ■V)Ey^{x - qp) 

^ /3=1 

+w„ A {{Vc, - i;;3) ■ V)5^^(x - qp) 
+ ^"oq(^^q)~M XI / d^x pa{x - qa)Va A By^{x - q/s) 

^ /3=1 

= : Ra,lit) + Ra,2it) + Ra,3it) + Ra^iif) 

with all time arguments at time t. Firstly, 



\Ki{t)\ 



(5.49) 



for a = 1, . . . , and t G [0, Te'^/^] by ( |CT1 ) and ( |CT| ). Since 5^ (a;) = -i; A V0„(x), by QCTI ), 
78D, and (g) also 

|i?.,4(t)| <Ce9/'. (5.50) 
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What concerns Ra,2{t)^ we may repeat the calculation in ( p.36| ) to obtain 
W^^^p(t) := j (fxp{x-q^{t))WU^^(^t){x-qp{t)) 

= ^ j j d^xd^yp{x - qa{t))p{y - qi3{t))Vl^Cvp(^t){x - y), 
with Cv{x) from ( ^.STp . Since supjg[o^Te-3/2] \'^lvCvf3{t){x)\ < C(l + we get as before that 

and hence by (|5.41|) 



So finally we have to bound Ra^sit), and this relies on a similar argument. Here we have 



and sup^gjQ | V^Ci',3(t)(^) I ^ C{1 + This in turn yields 

|V^„;3(t)| < Ce^ t G [0,T£-=^/2], a ^(3. 
Using the explicit form of Ey{x) and By{x), as in (|;39D, we then get for t G [0, Te-^/S] 

\RaAt)\<Ce'\vo.{t)-vp{t)\<Ce'/', 
by ( CT) . Summarizing (|^), (|530|) , (|53l|) , and (|531 it follows that 

\Ra{t)\ < Ce'^/\ a = l,...,N, t G [0,Te-3/2]. 
Consequently, by (|^, ( ^1481) , and (|3|) for a = 1, . . . , and t G [0, T^-^/Z] 



i^aWI < C(e^+\M^it)\ + \R^{t)\) < C 



e^^"^ + ( sup max \vg 



max Cfl 

</3<Ar ' '"^ 



Choosing maxi<^<jv le^gl < e with e sufficiently small we hence obtain 

1 N. 



sup \ia{t)\<Ce'"^, 

te[0,Te-3/2] 



(5.51) 



(5.52) 



(5.53) 



This completes the proof of Lemma 2.1 



□ 



6 Appendix B: Proof of Lemma |St2 



Here we give the proof of Lemma p.2| . We verify e.g. (b). To compare the left-hand side to the right- 
hand side of the assertion, we will insert some additional terms and estimate the corresponding 
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differences Dj{t), j = 1,2,3, for t G [to,Te ^/^], where to = 4(i?<^ + C*e ^). First we introduce 

_ifc.[,,,_i,2,^]l sin|fc|r^ 



Di{t) = ij^dr I d'k\ip{k)\'e 



\k\ 



-k 



V^y^ dr / d:'k\'p{k:)\e 



-ik-[TVp—^T^V0] 



sin \k\T 



= -^^jj d''xd'yip{x)ip{y) 

X J dr + ipit - ^)] - b - 9/3(^)]) - '^r{[x - - [?/ - Tt^/j])}, 

(6.1) 

as follows through application of the Fourier transform, with = Qait) — Qisit), and iprix) = 
(47r|a;|)~^ for |a;| = r whereas iprix) = otherwise. We claim that for x,y & IR^ with |a;|, \y\ < 
and t G [to,Te~^^'^] there exists a unique tq = To{x,y,t,$^a/3) G [0,to] C [0,t] such that 



T"0 



[-^a/? + X-qf3(t- To)] -[y- qf3{t)] 



(6.2) 



To see this, observe with 0(r) = r-|[^„^+x-g/3(t-r)]-[?/-g^(t)]| that > ^(0) > -{2R^+C*e-^) 
and 6''(r) > 1 — C^\fe by (p.6|) and ( p.8| ). For e so small that 1 — > 1/2 we hence obtain 

0{to) > -{2R^ + C*e-^) +to/2 = 2C*£"^ This shows 9 has a unique zero tq G [0,i(;o]- Moreover 
( |6.2| ) together with ( |2.6| ) implies 



^0 > IWI - k - Qfsit - To)] -[y- qi3it)]\ > C^e'^ - 2R^ - C^v^tq, 
whence also tq > Ce^'^ for e small. Similarly, we find a unique ti = Ti{x,y, 1,^^/3) satisfying 

1 



Tl 



(6.3) 



with Tl having the same properties as tq. By definition of ipr we therefore may simply write 

D,(t) = ~ I I dVy^ix)^{y) vJto' - t^'). (6.4) 



To estimate this, we calculate from ( |6.2| ) that 

V^Tq-^ = -ro^[{[^af3 + x-qf,{t-To)]-[y-q(,it)] 

+ + x- q^it - To)] -[y- g/3(t)]) • vpit - ro)V5ro|, 

with an analogous expression for V^rf ^. Therefore 

^i{r,'-r{')\ < c{To\p{t-To)-qp{t-Ti)\[l + \vp{t-Ti)\\V^Ti\^ 

+ - ^r'l |[W + Qpit - n)] -[y- qt3it)]\[l + \vf,it - rOllV^ril 
+ro^\Mt-ro)-vp{t-Ti)\\V^Ti\ 

+ro->^(t-ro)||V5(ro-n)|). (6.5) 
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^From (|6.2|), ( |6.3| ), and according to the Taylor expansion 

qp{t -t) = qp{t) - TVf, + ^r^i;^ + 0{e'l\''), 



cf. Lemma it follows that 

< Cy^\ro - Til + Ce^To + n)\To - nl + 

whence 

recall Ce:^^ < tq,ti < to = Differentiating ( p^ ) and ( |6.3| ) w.r. to ^ = we moreover 

get IV^TqI + iV^Til = and after a longer calculation which we omit also |V5(ro — ti)\ < 

C(£3/2 + y^|V^(ro-ri)|),thus 

\V^{To~n)\<Ce'/'. 

Utilizing these estimates and Lemma ^TTl in ( |6.5|) , we consequently obtain \'V ^{t^'^ — ri^)\ < Ce^^"^. 
Hence (|6.4| ) yields 

sup \Di{t)\<Ce'^'^ (6.6) 

as desired. Next, with 



D2{t) = I \ dr f d^k\0{k)fe-'''-^"'' 



it may be shown in a a similar way that 



2\\ sin|A;|r^ 



\k\ 



sup \D2{t)\<Ce'^/^. (6.7) 

te[to,Te-3/2] 



Finally we need to compare Jq dT(. . .) to the infinite rfr-integral and thus let 

\k\ 



D,{t)=ij^dr I d'k\m\'e-''-'''^{l-tk-[TV,-lT%]-^-r'ik.v,f) '^^k. 

sin \k\T 



With the notation 



\k\ 

this may be rewritten as 

D,{t) = J d'k mk)f^- Vgifo - {v/3 ■ V^)V^K, + ^{vf, ■ Vi:)V^K2 - ^{vf, ■ V^fV^K^^ 
Thus we only need to estimate 



d'k\0ik)\'Kp = j d'k\^{k)\'e-'^<-^ r dr^-^^^T^ 

\k\ 



dr i)r{iap + X - y) tP , (6.8) 
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the latter equality follows analogously to (|6.1|). However, for \x\, \y\ < and t G [to,Te ^/^] we 
obtain in case r = \^ai3 + x — y\ from ( p. 61) the contradiction 

4{R^ + C*e-^) =to<t<T<2R^ + < 2R^ + C*e-\ 



This shows the term in (|6.8| ) is identically zero for t G [to,Te~^/'^], and thus D^it) = for t G 



[to,Te ^/^]. Together with (O) and ( |6.7| ) this completes the proof of Lemma |3.2| (b). □ 
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